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Abstract 

We prove the existence of static, spherically symmetric solutions 
of the stellar dynamic Vlasov-Poisson and Vlasov-Einstein systems, 
which have the property that their spatial support is a finite, spheri- 
cally symmetric shell with a vacuum region at the center. 

1 Introduction 

Large stellar systems such as galaxies or globular clusters can be described 
by a density function / > on phase space. If collisions among the stars 
are neglected, / satisfies the so-called Vlasov or Liouville equation, which is 
then coupled to field equations for the gravitational interaction. Depending 
on whether one chooses a Newtonian or a general relativistic setting, the 
resulting nonlinear system of partial differential equations is the so-called 
Vlasov-Poisson or the Vlasov-Einstein system, respectively. In the present 
note we are interested in time independent spherically symmetric solutions 
of these systems. We call these solutions static, since due to the spherical 
symmetry their average velocity vanishes everywhere. The Vlasov-Poisson 
system then takes the following form: 

v-VJ-V x U-V v f = 0, (1.1) 
l(rV)' = 47rp (1.2) 
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where 

p(r)=p(x) = I f(x,v)dv. (1.3) 



Here x,vEJR 3 denote position and momentum, r— \x\, ' denotes the deriva- 
tive with respect to r, f = f(x,v) must be spherically symmetric, i. e., 
f(x,v) — f(Ax,Av) for every rotation A e SO (3), p(x)—p(r) denotes the spa- 
tial mass density of the ensemble and U(x) — U (r) is the induced gravitational 
potential. We assume that all particles in the ensemble have the same mass 
which — like all other physical constants — is set to unity. 

Under the corresponding assumptions the Vlasov-Einstein system takes 
the form 

^=.d x f-VIT^p'--d v f=o, (i.4) 

e~ 2A (2rA'-l) + l = 8nr 2 p, (1.5) 
e- 2X (2rp' + l)- 1 = 87cr 2 p, (1.6) 



p(r)=p(x)= / \ / l + v 2 f(x,v )dv , (1.7) 



where 



pW=pW = /(£_!!) 2 / (x,.)- 7 i= (1.8) 

denote the spatial density of mass-energy and radial pressure, respectively. 
If X = r(sin6'cos0,sin6'sin0,cos6') then the spacetime metric is given by 

rfs 2 = _ e 2, dj 2 + e 2X dr 2 + r 2 {M 2 + sin 2 

As boundary conditions we require asymptotic flatness, i. e., 

lim A(r)= lim u(r) = 0, (1.9) 

r—*oo r— >oo 

and a regular center, i. e., 

A(0) = 0. (1.10) 
For the Vlasov-Poisson system the corresponding boundary condition is 

limf/(r)=0. (1.11) 
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All solutions of the above systems known so far have the property that the 
support of p contains a ball about the center; the only steady states where the 
support does not equal such a ball are the axially symmetric ones obtained in 
[|]. The purpose of the present note is to construct solutions whose support 
is a finite, spherically symmetric shell so that they have a vacuum region at 
the center. Given the fact that the dynamical behaviour of both the Vlasov- 
Poisson and the Vlasov-Einstein systems is far from being understood, static 
solutions with new structural properties are of interest in themselves. How- 
ever, there is also a more specific motivation for the present investigation: 



In |TI] the gravitational collapse of spherically symmetric solutions of the 
Vlasov-Einstein system is investigated numerically. Static solutions provide 
useful test cases for the corresponding numerical scheme. Since the center 
of symmetry is particularly difficult to handle, it is important to have static 
solutions both with matter and with vacuum at the center in order to as- 



sess the performance of the numerical scheme. As shown in [|I(J symmetric 
solutions of the time dependent problem which vanish near the center re- 
main smooth. Thus, staying away from the center avoids analytic as well as 
numeric difficulties. 

The way to construct such steady states is now described. Since the 
system is time independent, the particle energy must be a conserved quantity, 
since it is spherically symmetric, the same is true for the modulus of angular 
momentum. Indeed, the quantities 

E = E(x,v) = ^v 2 + U(x), L = L(x,v) = \xxv\ 2 (1.12) 

are constant along solutions of the characteristic equations 

x = v, v = — VU(x) 
of the nonrelativistic Vlasov equation ( |1.1| ), and 

E = E(x,v)=e>* {r) Vl + v 2 , L = L(x 1 v) = \xxv\ 2 (1.13) 



are constant along characteristics of the relativistic Vlasov equation (|L 
Therefore, the ansatz 

f(x,v) = $(E,L) (1.14) 

satisfies the corresponding Vlasov equation and reduces the system to the 
field equation(s), where the source terms p or p and p now become functionals 
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of U or [a, which are obtained by substituting the ansatz ( |1.14| ) into ( |1.3| ) 
or (|1.7| ), ( |1.8|) respectively. In passing we note that every static, spherically 
symmetric solution of the Vlasov-Poisson system must be of the form ( |1.14j ), 
cf. 0. For the Vlasov-Einstein system this result, usually referred to as 
Jeans' Theorem, is not established. 

One can easily see that p becomes a decreasing function of r if / is a 
function of the particle energy E only, the so-called isotropic case. Thus, to 
obtain a nontrivial solution with a vacuum region at the center, / must also 
depend on the angular momentum L, and it must vanish for L small, say for 
L<L for some L >0. Once a solution of the field equation(s) is obtained, 
which has a vacuum region at the center, the main difficulty is to show that 
the support of the solution is actually bounded and the solution leads to a 
model with finite mass 

M = y p{x)dx <oo; 

in the case of the Vlasov-Einstein system this quantity is the so-called ADM 
mass. Finiteness of mass and support are obtained as follows. We take an 
ansatz function $, depending on the parameter L in such a way that for 
L = known results give the existence of a solution with finite mass and 
finite support, in this case a ball about the center. Then a perturbation 
argument in L is used to show that these properties persist also for L > 
but small. The smallness assumption on Lq can then be removed by a scaling 
argument. The details of this procedure together with the precise statements 
of our results are given in the next section for the Vlasov-Poisson system, 
and in the last section for the Vlasov-Einstein system. 

Before we go into this, we give a brief overview of the literature on the 
Vlasov-Poisson and the Vlasov-Einstein systems, starting with the former. 
We restrict ourselves to the stellar dynamics case; the plasma physics case, 
where the sign in the Poisson equation is reversed, is omitted. Global ex- 
istence of classical solutions has been established in ||, cf. also ^ [Tj|. As 
far as the existence of stationary solutions of the Vlasov-Poisson system is 
concerned we mention pi ffl. The main result on the initial value prob- 
lem for the Vlasov-Einstein system is a global existence theorem for small, 
spherically symmetric data ||. Spherically symmetric steady states for the 
Vlasov-Einstein system are constructed in [|, [|. 
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2 The nonrelativistic case 



Throughout this section we fix two parameters fc,/GlR with 

1 7 

fc>-l, 1>-1, k + l + ->0, k<3l + ~. 

2 2 

We make the ansatz 

f(x,v) = c Q (E -E) k + (L-L ) l + (2.1) 



where .E and L are defined as in ( |1.12|) , (•)+ denotes the positive part of the 



argument, and c >0, E <0, L >0 . It is a straight forward computation 
to show that with this ansatz 

P(r)=r 2l g(u(r) + ^J, (2.2) 



where 
and 



g{u) :=CoC kl (EQ-u) k + l+2 



c M :=2 /+ §7r l 1 ^=ds C s l+ ^{l- sfds, 
Jo y\ — s Jo 



and we have to solve 

I(rV)' = 4vrr 2 ^^ + ||), r>0. (2.3) 

The exponents k and I are kept fixed, while the parameters cq, E , and L 
may vary during our argument. The following theorem is the main result of 
the present section: 

Theorem 1 Let M>0 and Rq>0. Then there exists a static, spherically 
symmetric solution (f,p,U) of the Vlasov-Poisson system (\1.J\) , (\L-2j , 
where f and p depend on U via and U G C 2 ([0,oo[) flC 2 (lR 3 ) is 

a solution of satisfying the boundary condition ( \1.U\ ), p G C 1 ([0,oo[) R 



C 1 (1R 3 ) has total mass M, and suppp= [Ri,Ro] for some Ri & [0,i2o[, where 
Ri>0 provided Lq>0. Instead of prescribing M>0 and Rq>0 one may 
also prescribe M>0 and Ri>0. If Lq = and 0^/<l/2 then the asserted 
regularity holds only on IR 3 \{0}. 
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Note that we identify spherically symmetric functions of x with the corre- 
sponding functions of r— \x\. 

Proof: Let us fix some E and Cq = 1, and consider L = first. Then for 
Uq(0) <E prescribed there exists a unique solution U of 



: (r 2 U'y = 4irr 2l g(U) 



cf. |T]; this solution need not satisfy the boundary condition ( |1 . 1 1 ) , but we 
will take care of that later. By jj]], Uq induces a steady state with finite mass 
and finite support, which means that for some i?o > we have Uo(r) > Eq for 
all r > Rq. For L > we define 

U Lo (r) = U Q (0), 0<r<r Lo 



where 



Ln 



r L 



I2(E -U (0)Y 

note that Ul (r) + j% > E on ]0,r io [, so the right hand side of the Poisson 
equation vanishes on that interval. Now extend this towards the right by 
the solution of ( |2.3| ) with UL Q {r Lo ) = U (0), U' Lo (r Lo ) =0. The latter exists on 
[0,oo [, which can be shown using |l|], but this also follows from the arguments 
below. Upon integrating the Poisson equation we find 



4-7T r T 

K(r) = - s 2+2l g(U (s))ds } 
r z Jo 



where the latter integral is zero for r<ri . Now observe that Uo(r) >Uq(0) 
and Ul (t) > Ua(Q) for r>0, and on the set [{7o(0),oo[ the function g is 
bounded and Lipschitz; recall that we assume k + l + l/2>0. For 0<r<r^ o 
we have 

C 



U' {r)-U' Lo {r) 



< — 



Jo 



f+2l ds = Cr 2l+l 



and for r > r^ 



U'o{r)-U' Lo {r) 



C r 
<Uo(r Lo ) + - 2 

r Jr Lr) 



-,2+21 



U (s) 
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+ % f s 2+2l \Uo(s)-U Lo (s)\ds 

r Jr LQ 

<CLy l ^ + - 2 f r s 2 + 2l \U (s)-U L() (s)\ds, 

r Jr LQ 

where e>0 is such that 1 + / — e>0; constants denoted by C may depend on 
Uq and Rq, but not on r or L , and may change from line to line. Thus 

\U (r)-U Lo (r)\ <CL* + C f \ f a 2+21 \U (a) -U Lo (a)\ dads 

Jr L(j S Jr LQ 

<CL £ Q + C f s 2l+1 sup \U (a)-U Lo (a)\ds, 

JO 0<a<s 

and the latter inequality holds for 0<r<R . By Gronwall's Lemma, 

|C/ (r)-C/L (r)|<C^, O^r^; 

note that s 2l+1 is integrable over this interval. In particular, Ul must exist 
at least on this interval for L small. We conclude that 

U Lo (R )>E 

for L Q > sufficiently small, and by monotonicity, 

U Lo (r) + ^>E , r>R . 

We fix a small L > and let Ri = r^ and U — Ul etc. Then U leads to 
a steady state with finite radius and finite mass. The regularity assertions 
are obvious in the case L >0. In the case L = the exponent / has to be 
restricted in such a way that p'(0) = U'(0) = and U"(Q) exists. 

Let V(r)=U(r) + j$[. Then V(Ri)=E and V'{Ri)<0, whence V(r) < 
Eq in a right neighborhood of Ri. Thus the induced mass density is nontrivial 
and Ri is the radius of the inner boundary of its support. If we define 
M(r) =4n jQS 2 p(s)ds then M(r) is increasing with M(r) > for r > R iy and 
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Thus there can be at most one value of r > Ri where V'(r) changes sign, and 
if we define Ro = inf {r > Ri\V{r) = E }, this set being nonempty by what we 
showed above, then R >Ri, and suppp = [R^Rq], i. e., the support of the 
steady state consists of a single shell and not several nested ones. 

So far the boundary condition ( |1 . 1 1|) need not be satisfied, but U(oo) = 
lim r ^ 00 ?7(r) > E exists, and by slightly abusing notation we can redefine 

U = U-U(oo), E = E -U(oo). 

This leaves the distribution function / unchanged, and in addition is 
now satisfied. 

Finally we note that if / is a solution of the static Vlasov-Poisson system 
with mass M and suppp= [Ri,Ro] then the function 

fx,^x,v)=-f 3 X- 1 f(-fx,-f\- 1 v) 

is also a solution, with mass M(X,j) = A 2 7~ 3 M and support of the spatial 
density equal to [Ri/^jRo/j]. Choosing 7 and A appropriately any prescribed 
value for R (or Ri) and M can be obtained. Obviously, the constants Co, E , 
and Lq in the original ansatz (|2.1| ) are changed by this scaling, but the 
boundary condition ( |1 . 1 1| ) is not. 



3 The relativistic case 

Throughout this section we fix two parameters fc,/GlR with 

1 7 

/c>0, l> — , k<3l + -. 
~ ' 2' 2 

We again make the ansatz 

f(x,v) = c (E -E) k + (L-L o y + (3.1) 

where E and L are now defined as in (|1.13j ), and Co, -E'o>0, L >0. With 
this ansatz 



p(r) =r x e -W+^ g e^Jl + Lo/r 2 , (3.2) 



p(r) = r 2l e- {2l+A) »h e^Jl + L /r^ , (3.3) 
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where 



POO 

g(u):=c oCl {E Q -E) k + E 2 {E 2 -u 2 ) l+l ' 2 dE ) (3.4) 

J u 

Ku) :=^JjE -E) k + (E 2 -u 2 y^ 2 dE, (3.5) 



and 



Q:=27r/ -j==ds. 



Taking into account the boundary condition ( l.lOj) we can integrate the field 
equation (|1.5|) to obtain 



r Jo 



and substituting this into (|1.6|) reduces the static, spherically symmetric 
Vlasov-Einstein system to the equation 

/ 87]- rr \ — 1 / 4yj- pr \ 

p'(r) = fl -— J s 2 p(s)ds) UiTrp(r) + — J s 2 p(s)ds), (3.6) 

where p and p are now functionals of p given by (|3.2|), ( [3.3|) , (|3.4j ), (|3.5|) . The 
following theorem is the main result of the present section: 

Theorem 2 There exists a static, spherically symmetric solution 
(f,p,p,X,p) of the Vlasov-Einstein system j\1.4 ), ( \l-3j ), 
where f, p, and p depend on p via (p*. (\3.2( ), and p. j[ ) in a neighbor- 
hood of their support. \,p e C 2 ([0,oo[) nC 2 (IR 3 ) satisfy the boundary condi- 
tions firm, §rn\ ), and p is a solution of flFQ). p,p<E C^QO, 00 [) nC^K 3 ) 
with suppp = suppp= [Ri,Ro] for some 0<Ri<Ro<oo, where Ri>0 pro- 
vided Lq>0. The ADM mass M is finite, and one can prescribe M>0 or 
Ro > or Ri>0. If L = and 7^ / < 1/2 then the asserted regularity holds 
only on IR 3 \{0}. 

Proof. Consider first the case Lq = 0. As was shown in || there exists Eq > 
and a solution p of (|3l| with e w(0) <E and e w(i?o) > E for some Rq > 0. 
We choose R Q such that E Q < e^ ^ ) < E + 1. For any L > we define 



^2 e -2 M (o)_ 1 - 
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Then 



l + L /r 2 =E , e^Jl + L /r 2 >E , re[0,r Lo [, 



,W)(0) 



which means that PL (r) = po(0) solves (|3.6| ) on [0,r£, ] with pL ( r ) =PL (r) = 
0; in what follows pl and pi are always given in terms of pl by ( |3.2[ ) and 
( |3.3|) respectively By || Thm. 3.1] pi can be extended as a solution of 
( |3.6| ) for r>r Lo . We want to show that exp(p Lo (R )) > E for L >0 small 
so we may assume that pl (Ro) < Po(Ro) + 1 since otherwise we are done. By 
monotonicity, 

M0)<MO, p Lo (r)<p (R ) + l, re[0,R }. 

The functions g and h can be shown to be continuously differentiate, cf. || 
Lemma 2.1], and they vanish for u>E . Thus 



\pL (r) - p (r)\ < Cr 



21 



-(2Z+4) Mi0 



-(2^+4)^0 



g (e^oJl + L /r< 



g[e^oJl + L /ri\-g(e 



< C (r 2 Vo(r) -p (r)\ + r 2l - 1 ^j , r e]0,R 



and similarly 



(3.7) 



\PL (r)-p (r)\ <C [r \p Lo (r) - Po(r)\ +r 



L , relO,^ 



(3.8) 



Constants denoted by C may depend on /i and R , but never on r or L , 
and may change from line to line. From ( |3.6| ) we obtain the estimate 



l/4>)-/4Mi< 



1 8?r r 2 mj 
~ Jo s pL °( s ) 



+ 4vr 1- 



57T 



r Jo 



47rrp Lo (r) + 

s 2 p (s)(is 



1 87r f 2 r \ j 
1 / s p {s)ds 

r Jo 

s 2 p Lo (s)ds 



Air 
r 2 Jo 



r PL (r) + — s 2 p Lo (s)ds-rp (r)-— s 2 p (s)ds 
r z Jo r z Jo 
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Suppose that 

sup \fi Lo {r) - fi {r)\ <7 

0<r<_R 

for some 7>0 and Re]0,Ro]. Then by Q, 

y jT * 2 K (a) - Po (*) I ds < ° jT (>+ 2 7 + S 2 < +1 ^) da < Ci (7 + 

recall that / > — 1/2. Now we choose 7 >0 such that 

87r r 

0<r<Ho r "'O 

For Lq < 7 2 this implies that 



sup — / s 2 p (s)ds + 2Ci'y <1. 
i^r<Rn r JO 



0<r<i?, 

Thus 



8tt r r 87r /" r 

sup — / s 2 pL (s)ds < sup — / s 2 po(s)ds + 2Ci'j <1. 

<r<R T Jo 0<r<R T Jo 



1-— J s 2 p Lo (s)dsJ <C,0<r<R, 
provided (|3.9|) holds. This allows us to estimate the term I above: 

I<Cr 21 [ r s i\p Lo (s)-p (s)\ds 
Jo 

< Cr Al+2 JT + Cr 21 /> +2 |/i io ( S )-// (s)|d S . 
v Jo 

As to the second term we find that 

< Cr\p Lo (r) - Po ( r )\ + - f s 2 \p L (s) - p (s)\ds 

r z Jo 

<C v ^r 2 ' + C'r 2i+ VLoW-^o(r)| + C'|^ 2 VL (s)-/io(s)|^. 

Thus 

\ r i Lo (r)-p; (r)\<c(JL~ + J o s 2l \fi Lo (s) - fi (s)\ds 

on [0,i?], and by Gronwall's lemma, 

\fi Lo {r)-pio(r)\<cJT , re[0,R]. 
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By choosing Lq small we can make sure that ( |3.9| ) holds on [0,i?o] so that the 
previous estimate holds on [0,Ro\ as well. In particular, exp(p,L (Ro)) > E 
provided L is sufficiently small. We fix a sufficiently small L and let Ri = ri 
and /i = i^l etc. 

It is easy to see that \Jl + L /r 2 < E and thus p(r)>0 on some 
interval ]Ri,R[. Take Ro>Ri the smallest such R with the property that 
p = in a right neighbourhood of i?o- It is not clear that 1 + Lq / r 2 > Eq 

for all r > _R , but we can simply extend the solution towards the right of Ro 
by the corresponding vacuum solution. Thus, while /, p, and p depend on fi 
via (|3.1j) , (|3.2|) , and ( |3.3|) in a neighborhood of their support, this need not 
be true for all values of r. Clearly, /x(oo) = hm r __ >00 //(r) exists. If we redefine 

n(r)=fA(r)-n((x>), E = E e~^ 

we satisfy the boundary condition at infinity while the constant Cq in ( |3 . 1| ) 
is multiplied by e kfJ "(°°\ 

Finally, if f(x,v) defines a steady state, so does 

fa(x,v) = a 2 f(ax,v) 

for any a>0. The rescaled function f a has spatial support [a _1 _Rj,a _1 _Ro] 
and ADM mass 

a _i y f(x,v)\ / l + v 2 dvdx 

which shows that by rescaling a given solution we can get any prescribed 
value for the ADM mass, or the inner, or the outer radius. 
Final remark: It would be desirable to have a complete parametrization of 
all steady states as constructed in Thms. 1 and 2 (for k and I fixed), that 
is to say, a result of the form: For every 0<Ri<Rq<oo and M > there 
exists a unique steady state of the form fl2.1| ) or ( |3.1D with support [Ri,Ro] 
and mass M. Such a result can be obtained in the Vlasov-Poisson case for 
£o = 0, cf. [§. 
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